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Abstract 

In this Note we present the main results from the recent work [15], which answers several conjectures raised 
fifty years ago by Kac [9]. There Kac introduced a many-particle stochastic process (now denoted as Kac's master 
equation) which, for chaotic data, converges to the spatially homogeneous Boltzmann equation. We answer the three 
following questions raised in [9]: (1) prove the propagation of chaos for realistic microscopic interactions (i.e. in our 
results: hard spheres and true Maxwell molecules); (2) relate the time scales of relaxation of the stochastic process 
and of the limit equation by obtaining rates independent of the number of particles; (3) prove the convergence of the 
many-particle entropy towards the Boltzmann entropy of the solution to the limit equation (microscopic justification 
of the _H"-theorem of Boltzmann in this context) . These results crucially rely on a new theory of quantitative uniform 
in time estimates of propagation of chaos. 

Resume 

A propos du Programme de Kac en Theorie Cinetique. Dans cette Note, nous presentons les resultats 
principaux du travail recent [15], qui repond a plusieurs conjectures proposees il y a une cinquantaine d'annees par 
Kac [9]. Dans ce travail Kac introduit un processus stochastique a grand nombre de particules (aujourd'hui appele 
equation maitresse de Kac) qui converge, pour des donnees chaotiques, vers I'equation de Boltzmann spatialement 
homogene. Nous repondons aux trois questions suivantes soulevees dans cet article : (1) prouver la propagation du 
chaos pour des processus de collision realistes (dans notre cas : spheres dures et « vraies » molecules maxwelliennes) , 
(2) connecter les vitesses de relaxation du processus stochastique et de I'equation limite en obtenant des taux 
independants du nombre de particules, (3) prouver la convergence de I'entropie en grand nombre de particules vers 
I'entropie de Boltzmann pour la solution de I'equation limite (justification microscopique du theoreme H dans ce 
contexte). Tons ces resultats font appel de maniere cruciale a une nouvelle theorie d'estimations quantitatives et 
uniformes en temps de propagation du chaos. 



Version frangaise abregee 

Le programme de Kac en theorie cinetique consiste a comprendre comment deduire I'equation de Boltz- 
mann spatialement homogene a partir d'un processus stochastique de saut sur I'espace des vitesses a grand 
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nombre de particules. Le but de ce programme est de comprendre la notion de « chaos moleculaire » dans 
un cadre plus simple que celui de la dynamique complete des particules, ainsi que de donner une justifica- 
tion microscopique au theoreme H (croissance de I'entropie) et au processus de retour vers I'cquilibre. Nous 
renvoyons a la version complete pour I'introduction de I'equation de Boltzmann (1), des processus de saut 
consideres, ainsi que pour les definitions de la notion de chaos et de la distance de Wasserstein W\. 

Theoreme 1 (Resume des resultats principaux) On considere d > 2 et une distribution initiale fo G 
P{M.'^) n L°° a support compact ou possedant suffisamment de moments polynomiaux homes, et que Von 
suppose centree sans perte de generalite. Soit ft la solution correspondante de I'equation de Boltzmann (1) 
pour les spheres dures ou les molecules maxwelliennes (sans troncature angulaire), et soit la solution du 
processus de saut d N particules correspondant, avec pour donnee initiale : soit (a) la tensorisation de 
ff^ de fo, ou (h) la tensorisee f^'^ conditionnee a la sphere (definie par (2)). 

(i) Propagation de chaos quantifiee et uniforme en temps ; On considere le cas (a) ou fo = /®^. 
Alors 



avec a(N) — > lorsque N — > oo, et oil Heg^ designe la £-marginale d'une probabilite g sur (R'')"'^. 

(ii) Propagation du chaos entropique ; On considere le cas (b) ou f^ est conditionnee a, . Alors 
la solution est entropie-chaotique : 

(voir (3) pour les definitions des fonctionnelles H) avec 7 la probabilite gaussienne centree d'energie 
£ egale a I'energie de fo et la mesure de probabilite uniforme sur <S^. Cela fournit une derivation 

microscopique du theoreme H dans ce contexte. 

(Hi) Taux de relaxation independants du nombre de particules ; On considere le cas (b) oil f^ est 

conditionnee a . Alors 

W in. Yl i'y^W 

yN>l,yi<£<N,yt>0, ^ <p{t) avec p{t) ^ o, t ^ o. 

Dans le cas des molecules maxwelliennes, et si la donnee initiale possede une information de Fisher 
finie (voir (4)), on prouve egalementM N > 1, < jf H {ff\j^) < P{t) avec I3{t) -J> 0, i -> 0. 



1. Introduction 



Motivated by the understanding of irreversibility and "molecular chaos" in the Boltzmann equation [13,1,2] 
Mark Kac proposed in 1956 [9,10] the simpler and seemingly more tractable question of deriving the spatially 
homogeneous Boltzmann equation from a many-particle jump process, and he introduced a rigorous notion 
of molecular chaos in this context. He proposed the first proof of the propagation of chaos along time for a 
simplified collision process for which series expansions of the solution are available, and he showed how the 
many-particle limit rigorously follows from it. A key motivation for Kac was the microscopic derivation of the 
i?-theorem (monotonicity of the Boltzmann entropy) in this context which has remained open so far. Kac 
also raised the natural question of connecting the asymptotic behavior of the many-particle process and that 
of the limit nonlinear equation. In his mind this program was to be achieved by understanding dissipativity 
at the level of the linear many-particle jump process and he insisted on the importance of estimating its 
rate of relaxation. This has motivated beautiful works on this "Kac's spectral gap problem" [3,12,4], but 
so far this strategy has proved unsuccessful in obtaining relaxation rates which do not degenerate in the 
many-particle limit, see the interesting discussion in [5]. 

In this Note we present the main results in [15]. In this paper we develop a quantitative theory of mean- 
field limit which strongly relies on detailed knowledge of the limit nonlinear equation, rather than on detailed 
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properties of the many-particle Markov process. As the main outcome of this theory we prove uniform in 

time quantitative propagation of chaos as well as propagation of entropic chaos, and we prove relaxation 
rates independent of the number of particles (measured in Wasserstein distance and relative entropy) . All 
this is done for the two important realistic and achetypal models of collision, namely hard spheres and true 
(without cutoff) Maxwell molecules. This provides a first complete answer to the questions raised by Kac, 
however our answer is an "inverse" answer in the sense that our methodology is "top-down" from the limit 
equation to the many-particle system rather than "bottom-up" as was proposed by Kac. 



1.1. The Boltzmann equation 



The spatially homogeneous Boltzmann equation reads 

^{t,v)^Q(fJ){t,v), v&W, t>0, (1) 



where cZ > 2 is the dimension and Q is defined by 

Q{9,f){v) = ^ j B{\v -v^\,cose){gif' + g'fi- g^f - gf*) dv^da, 

where we have used the shorthands / = f{v), f = f{v'), g^ = g{v^) and g'^, = g{v'^). Moreover, v' and v'^, 
are parametrized by 

2 2 2 2 

Finally, 6 G [0,7r] is the deviation angle between v' — v'^ and w — defined by cos^ = a ■ u, u = v — v^, 
u = and B is the Boltzmann collision kernel determined by physics (related to the cross-section 

T^iv — v*,cr) by the formula B = \v — Vt:\^)- 

Boltzmann's collision operator has the fundamental properties of conserving mass, momentum and energy 



d_ 
di 



f <l){v) dv = j Qif, f) ct>{v) dv = 0, <p{v) = 1, V, \v\\ 



and satisfying the so-called Boltzmann's H theorem which writes (at the formal level) 

Rti R'' R'' 

where 7 is the gaussian with same mass, momentum and energy as /. Note that the H functional is the 
opposite of the physical entropy. 
We shall consider the foUlowing important physical cases for B (see [15] for more details) 

B = r{\v — v^\)b{cos6) with F, 6>0 given by one of the following formulas: 

(1) (HS) Hard Spheres collision kernel: B{\v — w*!, cos^) = F(|w — t;*|) = C\v — v*! for some C > 0. 

(2) (tMM) True Maxwell Molecules collision kernel: 

B(|w - I , cos 0) = 6(cos 0) -e^o C 0-^/^ for some C > 0. 

(3) (GMM) Grad's cutoff Maxwell Molecules kernel: B{\v - Uh.|,cos6') = 1. 



1.2. Kac 's program 



Kac 's jump process runs as follows: consider N particles with velocities fi , . . . , iiat € K"^. Compute random 
times for each pair of particles {vi,Vj) following an exponential law with parameter T{\vi — Vj\), take the 
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smallest, and perform a collision {vi,Vj) — )• {v*,Vj) given by a random choice of a direction parameter 
whose law is related to b{cos9), then recommence. This process can be considered on R''^, however it leaves 
invariant some submanifolds of Mf^^ (depending on the number of conserved quantities during collision) and 
can be restricted to them. In the original simplified model of Kac d = 1 (scalar velocities), the direction 
parameter is with collision rule 

V* = Vi cos^ + Vj sm6, Vj = —Vi sm.6 + Vj sva.6 

and the collision process can be restricted to Ei^~^{\/£N) the sphere with radius VSN, for any given value 
of the energy £. For the more realistic hard spheres of Maxwell molecules models, d = 3, the direction 
parameter is ct e §^ with collision rule 

Vi+Vj \Vi-Vj\ ^ Vi+Vj \Vi-Vj\ {Vi-Vj) 

V* = + — a, Va = — — (7 with a ■ = cos 6 

2 2 ' ^ 2 2 \vi-Vj\ 

and the collision process can be restricted to the sphere 

:-N ._^dN-ll ./^JC\ r. i.. , 1.^^=0}. (2) 

Kac formulated the notion of propagation of chaos that we shall now explain. Consider a sequence (/^)iv>i 
of probabilities on W^^ : the sequence is said f- chaotic if ^ f^^ when TV oo for some given one- 
particle probability / on W^. The meaning of this convergence is the following: convergence in the weak 
measure topology for any marginal depending on a finite number of variables. This is a low correlation 
assumption. It was clear since Boltzmann that in the case when the joint probability density of the 
iV-particle system is tensorized during some time interval into N copies f^^ of a 1-particle probability 
density, then the latter would satisfy the limit nonlinear Boltzmann equation during this time interval. In 
general interactions between particles prevent any possibility of propagation of the "tensorization" property, 
however if the weaker property of chaoticity can be propagated along time in the correct scaling limit it is 
sufficient for deriving the limit equation. Kac hence proved the propagation of chaos (with no rate) on the 
simplified collision rule above (with B = 1). His beautiful combinatorial argument is based on an infinite 
series "tree" representation of the solution according to the collision history of particles, and a Leibniz 
derivation-like formula for the iterated A/'-particle operator acting on tensor products. 

He then raises several questions that we schematize as follows: 

(i) The first one is concerned with the restriction of the; models as compared to realistic collision processes: 
can one prove propagation of chaos for the hard spheres collision process? 

(ii) Following closely the spirit of the previous question it seems to us very natural to ask whether one 
can prove propagation of chaos for the true Maxwell molecules collision process? This is 
related with long-range interactions and fractional derivative operators. 

(iii) Kac conjectures the propagation of the convergence of the A'^-particle il-functional towards 
the limit iJ-functional of the solution to the limit equation along time in the mean-field 
limit. Since the latter always decays for a many-particle jump process, in his words "// the above 
steps could he made rigorous we would have a thoroughly satisfactory justification of Boltzmann's H- 
theorem." 

(iv) He finally discusses the relaxation times, with the goal of deriving relaxation times of the limit equation 
from the many-particle system. This imposes to have estimates independent of the number of particles 
on this relaxation times: can one prove relaxation times independent of the number of particles 
in Wasserstein distance and/or relative entropy? 

This paper is concerned with solving the four questions outlined above. 
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2. Main results 



2.1. A few words on previous results 



For Boltzmann collision processes, Kac [9] has proved the propagation of chaos in the case of his baby 
one- dimensional model. It was generalized by McKean [14] to the Boltzmann collision operator for "Maxwell 
molecules with cutoff", i.e. the case (GMM) above (see also [19] for a partial result for non-cutoff Maxwell 
molecules) . Griinbaum [8] then proposed in a very compact and abstract paper another method for dealing 
with hard spheres, based on the Trotter-Kato formula for semigroups and a clever functional framework. 
Unfortunately this paper was incomplete for several reasons (see the discussion in [15]). A completely different 
approach was undertaken by Sznitman in the eighties [17,18] and he gave a full proof of propagation of chaos 
for hard spheres by a probabilistic (non-constructive) approach. Let us also emphasize several quantitative 
results on a finite time interval by Graham, Meleard and Fournier for Maxwell molecules models [11,6,7], 
and the works on the so-called "Kac's spectral gap problem" [3,12,4,5]. 



2.2. Main results 



Theorem 2.1 (Summetry of the main results) Consider some initial distribution fo G P(M.'^) n L°° 
with compact support or polynomial m,oment bounds, taken to be centered without loss of generality. Con- 
sider the corresponding solution ft to the spatially homogeneous Boltzmann equation for hard spheres of 
Maxwell molecules, and the solution fl^ of the corresponding Kac's jump process starting either (a) from 
the tensorization f®^ of fo or (b) the latter conditionned to (defined in (2)). 
The results in [15 J can be classified into three main statements: 

(i) Quantitative uniform in time propagation of chaos (with any number of marginals); 

VAr>l, Vl<^<iV, sn^^l^^^llA!Dl<a{N) 

f>0 " 

for some a{N) as N ^ oo, where HeQ^ stands for the i-marginal of an N -particle distribution 
, and where W\ is the Wasserstein distance between probabilities on W'-^: 

W\{p\,p2) '■= sup / (fi{dpi—dp2) where [j-Lip denotes the Lipschitz semi-norm. 

MLip(8«)<l^i 

In the case (a) f^ = f^^ one has moreover explicit power law rate (for Maxwell molecules) or 
logarithmic rate (for hard spheres) estimates on a. 

(a) Propagation of entropic chaos.- Consider the case (h) where the initial datum of the many-particle 
system is restricted to . Then if the initial datum is entropy- chaotic in the sense 

1 

TV 



with H{f,^\j^) ■■= j ^ \og^^^{dV) and H {f,\^) := j fo log^dv (3) 

and where 7 is the gaussian equilibrium, with energy £ and 7^ is the uniform probability measure on 
, then the solution is also entropy- chaotic for any later time: 

Vt>0, liJ(/f|7^) ^iJ(/,|7), TV ^+00. 

Since our /g^ is entropy- chaotic, this proves the derivation of the H-theorem in this context. 
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(Hi) Quantitative estimates on relaxation times, independent of the number of particles; Con- 
sider the case (b) where the initial datum of the many-particle system is restricted to . Then 

yN >i, yi<e< N, yt>o, ^^i^^ft ^^^eii )) ^^^^^ with p{t) ^ o, t ^ o. 

Moreover in the case of Maxwell molecules, and assuming moreover that the Fisher information of the 
initial datum fo is finite: 

[\l^dv<+oo, (4) 
J To 

the following estimate also holds: ViV > 1, < H (/^^It^) < l3{t) with l3{t) -J> 0, t ->■ 0. 

3. A few words on the methods and proofs 

Let us briefly explain some ideas underlying the result Theorem 2.1-(i). The other results are then obtained 
on the basis of this key estimate, combined with the other latest results obtained in this field. 

- We aim at reducing the problem to a stability analysis when approximating a linear semigroup. To this 
purpose a key idea is to compare the linear iV-particle dual evolution in Ci){R'^^) with the (linear!) 
push-forward evolution associated with the limit equation. 

- This push-forward semigroup is defined as follows: if S^^ denotes the nonlinear semigroup of (1), this 
push-forward semigroup is defined on P(P(M'')) by T~[$](/) = 

- In order to make this comparison between semigroups, we use the empirical measure fXy = (X^^i J^J/AT 
in order to embed the dynamics in P(]R''^) into a dynamics in P(P(M'')). 

- One then considers the following term to be estimated 

for some test function ip only depending on £ variables. 

- The approximation of these marginals by empirical measure estimates yields a first error term on the 
A/'-particle semigroup 

|(5f(/o^),<p0l«^-O-(5f(/o^),4o^{)f)| with Rf^if):= I ^f^\dv,...dv,) 
which is controlled by combinatorial arguments, and then a second error term 

which is controlled thanks a stability for measure solutions of the limit equation. 

- Finally there remains the most important term where the two dynamics are effectively compared 

This term is controlled by using (1) a quantitative argument a la Trotter-Kato in order to express the 

difference of semigroups in terms of the difference of their generators and G°°, (2) a consistency 
estimate between those generators, (3) a stability estimate on the limit equation. 

- There is a loss of derivative in the consistency estimate in the sense of differentiable functions acting 
on P(R''), which lead us to develop a differential calculus on this space adapted to our purpose. 

- The stability estimate means, once translated on the original nonlinear semigroup S^^ of (1), a propa- 
gation of a bound C^+^(P(]R'^)) on S^^. The role played by such stability estimates is a key novelty of 
our study. Proving them for Boltzmann is also one of the most technical aspects of [15]. 

- There are many possible choices of distances on the space of probabilities (total variation but also many 
non-equivalent weak measure distances), and it is a crucial point that our method is flexible enough to 
allow for many such different choices adapted to the equations it is applied to. 
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